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Abstract 



We describe the results that have so far been obtained in the classifica- 
tion problem for integrable (2+l)-dimensional systems of hydrodynamic 
type. The Gibbons-Tsarev (GT) systems are most fundamental here. A 
whole class of integrable (2+l)-dimensional models is related to each such 
system. We present the known GT systems related to algebraic curves of 
genus g = and g = 1 and also a new GT system corresponding to al- 
m ' gebraic curves of genus g = 2. We construct a wide class of integrable 

models generated by the simplest GT system, which was not considered 
previously because it is "trivial." 



Keywords: dispersionless integrable system, hydrodynamic reduction, sys- 
tem of Gibbons-Tsarev type 
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1 Introduction 

Integrable systems play an important role in both mathematics and physics. 
Unfortunately, a rigorous universal definition of integrability applicable to dif- 
ferential equations of all types is lacking. Different views of the integrability 
problem were described, for example, in [T]-[3]. 
Here, we consider systems of the form 

^ Q ^ Q ^ Q 

^ a ^ U )^ + ^ b ij( U )-Hr + ^2 C ij( U )-HZ= > U=( Ul ,...,U n ), (1) 

3 = 1 3 = 1 U 3 = 1 

where i = 1, . . . , n + k and k > 0, and also equations that are reducible to the 
above form by substitutions, for instance, quasilinear equations 

A\Zu + AiZyt + A$Z x t + A^Zyy + A§ Z xy + AqZ xx — 0, (2) 

where Ai = Ai(Z x ,Z y ,Zt), and equations of the form 

! Zyt, Z X t, Zyy, Z X y , Z ' XX ) 0. (3) 
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Here and hereafter, we use the notation Z x = dZ/dx, etc. To include equations 
and systems of types ([2]) and ([3]) in our consideration, we must allow the number 
of equations in system (fTJ) to exceed the number of unknowns, i.e., we must be 
able to consider the case k > 0. 

The existence of a dispersionless zero-curvature representation [5], [B] (also 
see [7]) or (which is the same) a pseudo-potential representation is a fundamental 
property of integrable systems (fTJ . This means that system (fTJ is equivalent to 
the compatibility conditions for the pair of the Hamilton- Jacobi equations of 
the form 

$ y = A($ x ,ui,...,u n ), $ t = B($ x ,ui, . . . ,u„), (4) 
where $(x, y,t) is a scalar function. As an example, we consider the system 

u y = v x , Vy = u t — uu x . (5) 

We can easily verify that it admits the pseudopotential representation 

<S>y = ^+U, $ t = -g. + U<S> X + V. (6) 

Eliminating the unknown v, we can write system ([5]) as the equation 

Uyy = (U t ~ UU X ) X , 

which is known as the dispersionless Kadomtsev-Petviashvili (KP) equation or 
as the Khokhlov-Zabolotskaya equation. 

For some integrable models, the functions A and B in the pseudopotential 
representation depend on an additional spectral parameter A. An example is 
the dispersionless Hirota equation 

CLiZ x Z yt + a 2 Z y Z xt + a 3 Z t Z xy = 0, a\ + a 2 + a 3 = 0, (7) 

which has the pseudopotential representation 



a 2 X Z v Z t 

> y = -7=r®x, $t = A — ■ 

a\X + 03 Z x Z x 



As in the case of integrable equations of the Korteveg-de Vries (KdV) equation 
type, the dependence on A can be polynomial, rational, and so on. Assuming the 
existence of a pseudopotential representation with a prescribed A-dependence, 
we can construct examples of integrable systems ((TJ). For small n, we can also 
classify systems admitting pseudopotential representations without the spectral 
parameter. But the very existence of a pseudopotential representation cannot 
serve as a universal integrability criterion for systems (TTJ because various types 
of pseudopotential representations exist. 

The integrability of system ([1]) is often related to its representation as the 
commutation condition for a pair of vector fields [8]-|ll|. which may or may 
not depend on the spectral parameter A. If a pseudopotential representation 
exists, then these vector fields are Hamiltonian. Then, apparently, we cannot 
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select a constructively described class of vector fields covering all the known 
examples and provide a rigorous universal definition of integrability based on 
the commutation of vector fields. 

In our view, we currently have only one property of integrable systems ([TJ 
that can be taken as a universal constructive integrability criterion, which means 
that we not only can verify this property for a given system but also can use it 
to construct new integrable models and to provide a complete classification of 
these systems in lower dimensions. This integrability criterion is that system |T]) 
has a sufficient number of the so-called hydrodynamic reductions. The notion 
of a hydrodynamic reduction arose in |12j . where the Benney chain reductions 
were found. It turned out that iV-component hydrodynamic reductions can 
be described in terms of the compatible overdetermined Gibbons-Tsarev (GT) 
system 

^ u = ^r^' i,j = l,...,N, (8) 

Here, the unknown functions <p\, . . . ,Pn, and u depend on some variables r 1 , . . . , r 
and di — d/dr l . The same system ((5J) appears when describing hydrodynamic 
reductions of the Khokhlov-Zabolotskaya equation. 

It was mentioned in the pioneering papers [13]— [T5] that the presence of hy- 
drodynamic reductions can be effectively used to classify integrable systems ([T]) . 
Below, we briefly list the main classification results. We describe the hydrody- 
namic reduction method itself in Sec. [2] 

The problem of classifying integrable systems of the form 

v t + av x + pv y + qw y = 0, w t + bw x + rv y + sw y = (9) 

was considered in |14j . A system of 16 nonlinear partial differential equations for 
the unknown coefficients a(v,w), b(v,w), p(v,w), q(v,w), r(v,w), and s(v,w), 
which was equivalent to the existence of hydrodynamic reductions, was written 
there. Because the obtained system for the coefficients is in involution, we 
can use it to study general properties of integrable systems ([9]). In particular, 
it was proved in |14) that any integrable system (|9]) admits a pseudopotential 
representation and that the space of conservation laws of hydrodynamic type is 
three-dimensional. The authors found many particular solutions of the system 
of equations for the coefficients and thus constructed new examples of integrable 
systems But the authors did not construct the general solution explicitly. 

We can easily verify that for n = I = 2, general system ([1]) is reducible to 
canonical form ([9]) by a coordinate change of the form 

v = ipi(ui, u 2 ), w = ip2(ui,u 2 ). (10) 

But fixing canonical form ([9]) destroys the group GL(3) that acts on the space 
of integrable systems ([1} by changing the independent coordinates: 

i = knt + k 12 x + ki 3 y, y = k 21 t + k 22 x + k 23 y, x = k 3 %t + k 32 x + k 33 y. (11) 
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Reducing the action of this group to its simplest form, Odesskii [TD] found 
another canonical form in which the system of equations for the coefficients can 
be solved explicitly in terms of generalized hypergeometric functions. 

A special subclass of integrable equations © with a Lagrangian of the form 
jC(Z t , Z x , Z y ) was considered in [IT], [T5]. Necessary and sufficient conditions 
for the Lagrangian integrability were obtained in |17j in the form of an overde- 
termined system of partial differential equations for the function C. General 
properties of integrable Lagrangians and new examples were also constructed. 
A new integrable Lagrangian C ~ Z t Z x Z y found there is the most interest- 
ing. General integrable Lagrangians were described in [IS] in terms of modular 
forms. 

The general classification problem for integrable equations ([2]) was consid- 
ered in |19) . where an overdetermined system in involution for the coefficients 
Ai(Z x , Z y , Z t ) (which is equivalent to the existence of hydrodynamic reductions) 
was written, general properties of integrable equations (J2J were investigated, and 
new examples of these equations were found. The general solution of this overde- 
termined system was found in [2D]. It depends on five arbitrary constants and 
can be expressed in terms of generalized hypergeometric functions [21] of three 
variables after an appropriate coordinate change. We note that in addition to 
depending on the five essential parameters, the set of integrable equations ([2]) 
depends on 15 more parameters generated by the action of the total linear group 
of transformations of the parameters Z, x, t, and y, which is the invariance group 
for equation class (0). 

Necessary and sufficient conditions for the integrability of equations of form ((3|) 
were obtained in [22]. It turned out that the function F does not contain es- 
sential parameters. Moreover, interesting classes of particular solutions of the 
system for the function F were found in [22]. The general solution was con- 
structed in [20] in terms of the generalized hypergeometric function. 

The classification problem for infinite integrable divergent chains of the form 

u u = Fi(ui,u 2 ) x , u 2t = F 2 (u ll u 2l u 3 ) Xl ...,u it = F i (u 1 ,u 2 , . . .,u i+ i) x , . . . 

(12) 

was considered in [23]-[25]. The integrability was understood in [23] as the 
existence of a commuting flow of the form 

uiy = Gi(ui,u 2 ,u 3 ) x , u 2y = G 2 {ui,u 2 ,u 3l U4) Xl . . . ,u iy = Gi(ui,u 2 , . . . ,u i+2 ) x , ■ ■ ■ 

(13) 

A necessary condition for the existence of hydrodynamic reductions (the vanish- 
ing of the Haantjes tensor) was used in [25] to produce the system of equations 
for the functions F% and F 2 that must be satisfied for the chain to be integrable. 

In the above classification problems, in addition to having hydrodynamic 
reductions, "general position" integrable models have pseudopotential represen- 
tation (j4]) without the spectral parameter. All the pseudopotentials arising here 
are related to some rational curves. Krichever constructed [26], [6] a class of 
pseudopotentials with arbitrarily large n related to algebraic curves of an ar- 
bitrary genus n. It can be verified that the set of Krichever pseudopotentials 



5 



corresponding to a rational curve does not include all the examples found when 
classifying system (|9|) for n = 2. 

The definition of the pseudopotential integrability and the complete descrip- 
tion of integrable pseudopotentials for n — 1 were given in terms of hydrody- 
namic reductions in [27] • The answer was formulated in terms of hypergeometric 
functions. 

Numerous examples related to rational curves lead to the problem of con- 
structing a class of pseudopotentials with arbitrary n that contains all these ex- 
amples (possibly as reductions). Such a pseudopotential class was constructed 
in [TO] , [20J • The result was expressed in terms of a special class of generalized 
hypergeometric functions of n variables and contained n+2 essential parame- 
ters. When all these parameters are integers, the result can be expressed in 
terms of elementary functions. The problem of explicitly describing reductions 
of these potentials remains open. The obtained results were generalized to the 
case of an elliptic curve in [55]. Generalized elliptic hypergeometric functions 
appear in the answer in [29] ; the class of elliptic hypergeometric functions there 
is apparently new. 

Here, our presentation is based on the most fundamental notion of the hydro- 
dynamic reduction method: the GT systems. Although such systems appeared 
previously [30] -[31], a rigorous definition of the class of such systems, their 
equivalences, etc. (see Sec. [3J, is apparently still lacking in the literature. We 
present the known GT systems related to algebraic curves of genus g = and 
g = 1 and also a new GT system related to algebraic curves of genus g = 2. 
A whole class of integrable (2-(-l)-dimensional models is related to each of the 
GT systems. In Sec. IU we describe the results in ^20 , [28] in terms of the GT 
system. In Sec. [5] we construct a new broad class of integrable models generated 
by the simplest GT system, which was not considered previously because it is 
"trivial." These systems have a pseudopotential representation that depends on 
the spectral parameter rationally. 

In Sec.[0J we share our experience acquired during several years of work with 
systems (JJ). We devote Appendices IA1 [Dl to conservation laws, pseudopotentials, 
multidimensional integrable systems ([1}, and integrable chains. Appendix [Al is 
important for understanding the relation between integrable systems ([1} and 
integrable models of the type of Eqs. ([2]) and (J3j) . We describe the relation 
between the GT systems and the pseudopotential representation without the 
spectral parameter in Appendix |B] which is important both theoretically and 
technically. Examples of multidimensional integrable systems of hydrodynamic 
type are in Appendix [C] where we also briefly discuss how to generalize the 
GT-system notion to the multidimensional case. We devote Appendix ID] to the 
problem of classifying arbitrary (not necessarily divergent) hydrodynamic-type 
chains with the unit shift. Moreover, we present examples of integrable multi- 
dimensional chains in Appendix [D] We list unsolved problems in Appendix [El 
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2 Hydrodynamic reduction method 



It is well known that integrable equations of the KdV equation type 

-K-s SO <»> 

have families of explicit solutions depending on arbitrary constants c\, . . . , cjv for 
any AT. All these finite-gap and solitonic solutions can be constructed using the 
so-called ODE reductions. We call a pair of compatible A^-component systems 
of ODEs 

rl = f{r\...,r N ), r\ = ^(r 1 , . . . ,r N ), i = l,...,N, (15) 

the Af-component ODE reduction of Eq. (fT4]l if we have a function ^(r 1 , . . . , r N ) 
such that u — U(r 1 (x 1 1), . . . , r N (x, t)) satisfies Eq. (|14[) for any solution r 1 (x, t), . 
of system (fT5|) . Clearly, a solution u depends on N arbitrary parameters, which 
are initial values for system (|15[) at the general position point. We can take the 
existence of special classes of the ODE reductions with an arbitrary N as an 
integrability criterion for Eq. (fT4|). For instance, we can assume that Eq. (| 14[) 
admits a sequence of differential constraints of the form 



d m u „ / du d r ' 



= G r . 



u. 



dx m \ ' dx ' ' dx m ~ 

with an arbitrary m. Clearly, we can rewrite Eq. (fT4|) and such a differential 
constraint as a pair of compatible dynamical systems in the variables x and 
t. Another example of an ODE reduction is the Dubrovin equations for zeros 
of the Baker-Akhiezer function in the KdV equation theory. But much more 
effective and constructive integrability criteria exist in the two-dimensional case. 
For example, we can take the existence of higher local commuting flows (higher 
infinitesimal symmetries) or the existence of higher conservation laws (see [35] 
and the references therein) as such a criterion. 

Higher local symmetries and/or local conservation laws are absent if the 
number of independent variables is d > 2 (see [36] for generalization of the 
symmetry approach to the case of nonlocal symmetries). The existence of TV- 
component reductions can be considered one of the most prospective approaches 
for seeking new integrable models in this situation. We note that instead of 
the ODE reductions, we must there take some compatible systems of partial 
differential equations of dimension not higher than d — 1 as reductions. 

Such an approach was applied to systems of form ([1]) in [13] . where pairs of 
compatible semi-Hamiltonian (see (|18j) below) systems of hydrodynamic type 

rj = A < (r 1 ,...,r Ar )r« ) rj = M V, • • • ,r*>*, i=l,...,N, (16) 

were taken as reductions. By the definition of the reduction, the correspond- 
ing solutions of system ([T]) are provided by some functions u % (r 1 , . . . , r N ) , i = 
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1, . . . , n, that transform any solution of system (TTH1) into solution (TTJ). Such so- 
lutions describe the interaction of N plane waves in hydrodynamics. They are 
sometimes called the TV-phase solutions. 

Clearly, the generalized solution of a compatible overdetermined system of 
form (I16p contains N arbitrary functions of one variable. It turns out that 
the functions A 1 and fi 1 in reduction (|16p may contain additional functions of 
one variable as functional parameters; the number of these functions does not 
exceed N. The existence of hydrodynamic reductions (fTS)) locally parameterized 
by N functions of one variable with arbitrary N was taken to be an integrability 
criterion for systems of form (fTJ) in [T3] . The corresponding TV-phase solutions 
depend on 2N arbitrary functions of one variable. 

The approach based on hydrodynamic reductions is universal. This means 
that all the currently known integrable systems (fTJ) admit hydrodynamic reduc- 
tions. All the notions in this approach can be rigorously defined (see below). It 
was demonstrated in [T3] that the existence of hydrodynamic reductions can be 
verified algorithmically and can be effectively used to classify integrable cases. 

We recall the definitions and basic notions of the hydrodynamic reduction 
method. Although they essentially coincide with those in [13], we introduce 
some technical improvements to make the definitions more rigorous. 

The geometric theory of integrable (l+l)-dimensional hydrodynamic-type 
systems of the form 

r\ = X{r\...,r N )ri, i = l,...,N, (17) 

was described in [37] . It was shown in [38] that these systems can be integrated 
using the generalized hodograph method. 

Definition 1 System (|17l) is said to be semi-Hamiltonian if we have the rela- 
tions 

9i^^- = di , 3 \ , di = ^- (18) 

\ l — A m A J — A m ' Qr l 

Semi-Hamiltonian systems have an infinite number of commuting flows and con- 
servation laws of hydrodynamic type [381, [39]. The main geometric object re- 
lated to semi-Hamiltonian system (jTTJ) is the diagonal metric gj.}., k = 1, . . . , N , 
where 

-ddogguu = — — A 
which is called the metric associated with (|17[) 



^logff fc fc = u l xfc , i^k, (19) 



Definition 2 A hydrodynamic reduction of system (TTJ is governed by a pair 
of compatible semi-Hamiltonian systems of hydrodynamic type (|16p and by the 
functions Ui(r 1 , . . . , r ), . . . , M n (r 1 , . . . , r ) such that the functions 

«i =u 1 {r\...,r N ), u n =u n {r 1 ,...,r N ) (20) 

satisfy system ([1]) for any solution of system (|16p . Following \13f -J15 f , we say 
that system (TTJ is integrable if it admits the maximum possible number of hydro- 
dynamic reductions. Namely, substituting (|20p in ([1}, eliminating derivatives 
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with respect to t and y using (|16p , and equating the coefficients of r x to zero, 
we obtain 

n n n 

^aij(u)X l diUj + ^bij(u)iJ, l diUj +'y]c i j(u)diUj = 0, (21) 

3=1 3=1 3=1 

i = 1, . . . , n + k, I = 1, . . . , N. 

For each I, we have a linear overdetermined system for n unknowns diui, . . . , diu. 
whose coefficients are independent of I. Because this system must have a non- 
trivial solution, all its nxn minors must vanish, which results in the system of 
algebraic equations for A' and fi , which is the same for all I. We assume that 
this system is equivalent to a single algebraic equation 

P(A',m') = (22) 

(otherwise, A' and u 1 are fixed, and we then have just a finite number of hydrody- 
namic reductions). Equation (|22[) determines the so-called algebraic dispersion 
curve. Let p be a coordinate on this curve. Then Eq. (|22p is equivalent to the 
equations 

A' = F(pi,Ui,...,Un), n l = G{pi,ui,...,u n ) 

for some functions F and G. It was assumed in [13] that linear system (|2Tj) has 
a unique solution up to proportionality at a general value of pi . Solving this 
system, we obtain 

diu m = g m (pi,ui, . . . ,u n )diUi, m = 2,...,n, i = l,...,N, (23) 

for some functions g rn . We rewrite (| 16|) in the form 

r\ = F(p i ,u 1 ,...,u n )rl, r l y = G(pi,ui, . . . ,u n )r l x , i = l,...,N. (24) 

It is easy to see that the compatibility conditions for (pM)) are 

FM-Ffa) G{ Pl )-G(p 3 ) 1 ] 

(where we omit the arguments u\, . . . , u n of F and G). We express diPj from (|25|) : 
dipj = f(pi,pj,ui,...,u n )diUi, i^j, i,j = l,...,N. (26) 
Finally, the compatibility conditions didjU m = djdiU m result in the equation 

d i djUi=h{pi,pj,u 1 ,.,.,u n )d i uidjUu i ^ j, i,j = 1, . . . ,N. (27) 
Collecting Eqs. (|2"3")l , (|2^|) . and (|2T|) together, we obtain the system of equations 
diPj = f(js>i,Pj,u 1 ,...,u n )diUx, i^j, i,j = l,...,N, 

diu m = g m (Pi,ui, . . . ,u n )diUi, m = 2,...,n, i = l,...,N, (28) 

d i djU 1 =h(pi,Pj,ui,...,u n )d i u 1 djUi, i ^ j, i,j = l,...,N. 
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Here, pi, ■ ■ ■ ,pn and m, . . . , u n are functions of r 1 , . . . , r , N > 3, and di — 
d/dr 1 . Because system ([T]) must have the maximum number of reductions by 
assumption, system (125)) must be in involution (i.e., it must be completely com- 
patible). This means that the compatibility conditions dkdjPj — djdkPj and 
dkdidjUi — didkdjiii are satisfied by virtue of system (|2"g|) . Solutions of sys- 
tem (|28p therefore depend on 2N functions of one variable. We can take these 
functions to be boundary values of the Goursat problem ui(0, . . . , 0, r l , 0, . . . , 0) 
and Pi(0, . . . , 0, r 4 , 0, . . . , 0), where i = 1, . . . ,N. 

3 Gibbons— Tsarev systems 

In the preceding section, we showed that a system of form (|28|) is related to 
each system ([1]) that has the maximum set of hydrodynamic reductions. 

Definition 3 A compatible system of form (|28[) is called an n-Eeld GT system. 

Definition 4 Two GT systems are said to be equivalent if they are related by 
a transformation of the form 

Pi -)• \(j)i,ui, . . . ,u n ), i = l,...,N, (29) 

u m ->• Mm •••.'"„), m=l,...,n. (30) 

Definition 5 Transformation (|29|> . p0[) is called a GT-system automorphism 
if it transforms this system into itself. 

It is easy to see that system (|28l) is compatible if and only if the functions 
/, gk, and h satisfy a system of functional equations independent of N. For 
example, for the one-field GT system 

diP] = f(pi,Pj,u)diU, did j u = h(p ll pj,u)d l ud :j u, i^j, i, j = 1, . • . , N, 

(31) 

these equations are 

fu{P2,P3,u) - f u (p 1 ,p 3 ,u) + fp 2 (p2,P3,u)f(p 1 ,p 2 ,u) ~ f pi (p 1 ,p 3 ,u)f(p 2 ,Pl,u) + 

+ fp 3 (P2,P3,u)f(pi,p 3 ,u) - f p3 (p 1 ,p 3 ,u)f(p 2 ,P3,u) + 
+ f(P2,P3,u)h(p 1 ,p 2 ,u) - f(p 1 ,p 3 ,u)h(p 1 ,p 2 ,u) 

h u (p2,P3,u) - h u (pi,p 3 ,u) + h P2 (p2,p 3 ,u)f(p 1 ,p 2 ,u) - h Pl (p!,p 3 ,u)f(p2,pi,u) + 

+ h P3 (p2,p 3 ,u)f(p 1 ,p 3 ,u) - h P3 (p 1 ,p 3 ,u)f(p2,p 3 ,u) + 

+ h(P2,P3,u)h(p 1 ,p 2 ,u) - h(p 1 ,p 3 ,u)h(p 1 ,p 2 ,u) 

(32) 

3.1. Examples of GT systems 



10 



Example 1 Let ao, a\, and 0,2 be arbitrary constants. Then the equations 

a 2 p 2 + axPj +a 
OiPj = diU 

I 1 - /'.. . . , Ar . . . 

t,j = 1,...,N, t^j, 

0,0,11- (pi-pjf ° lUa ^ 

determine a one-field GT system. The original GT system (see Sec. [7]) corre- 
sponds to the case where a 2 = &i = and ao = 1. In the general position case, 
we can use a linear transformation of the variables pi to reduce the system to 
the form 

d iPj = p -^i^ di u, 
Pi - pj 

i,j = l,...,N, i^j. (33) 

2piPj — pi — pj 

Example 2 Let P(x) — a 3 x 3 + a 2 x 2 + a\x + oq. Then the system of equations 
P{u)(Pi- P j) 

i,j = l,...,N, i^j, 

n n „ _ K 2 (pi,Pj)u + K^p^p^u + Kpjpupj) 
didjU ~ P(u)( Pi - Pj ) 2 % 3 1 

(34) 

where 



K 2 (pi,Pj) = 2a 3 (p 2 - pj) 2 , 

Ki(pi,pj) = ~a 3 (p 2 p 3 + p l p 2 ) + a 2 (p 2 +p 2 - ipiPj) - ai(pi +Pj)- 2a , 

Ko(Pi,Pj) = "2a 3 p 2 p 2 + a 2 (p 2 pj +P1P 2 ) + ai(pf + pf) + a Q (p l +p 3 ), 
is a one- field GT system. Using transformations of the form 

au + b api + b 

u ^ — n> P*~* — r^' 

cu + a cpi + a 

we can reduce the polynomial P to one of the three canonical forms P(x) = 
x(x — 1), P(x) = x, or P(x) = 1. We note that in the general position case 
where P{x) = x(x — 1), the system has an automorphism group S4 interchanging 
the points 0, 1, 00, and u and generated by the elements 

u Pi 

erf. u 1 - u, pi-tl — Pi, <j 2 : u -, pi -> 



u — 1 pi — a 

and 

(1 - u)pi 

cr 3 : u -> 1 - u, pi -> . 

Pi - u 
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Example 3 We set 

9(Z,T) = ^ ( _l)« e a 1 ri(«+«(a-l)r/ a ) i ^ r) = 

XTien i/ie system 

daP/3 = 77—(p(Pa-Pfs)~P{Pa))d a T, d a dj}T = p' ' (p a ~ pp)d a Tdf3T, (35) 

Z7TI TTl 

where a, f3 = 1, . . . , N, a ^ j3, is a one-field GT system. The field u here is 
the modular parameter r. We note that the same system completed with some 
equations with a — /3 appeared in in relation to a very different problem. 

Example 4 Let u, v, and w be the coordinates on the moduli space of curves 

y = x(x — l)(x — u)(x — v)(x — w) (36) 

of genus g = 2. Then the formulas 

v(v-l)(pi-u) 
diV = ~? 777 7 5 4 w, 

U(U ~ l)(Pi - V) 

u(u - l)(pi - w) 

g = PjiPj " iXgj ~ u)(Pi ~ v)(Pi -w) + y(pi)y(pj) g u 
3 u(u - \){jpi - v)(pi - w)(pi - Pj) 



didoU — 



«(Pt + Pj) + (u ~ l)(pf + Pj) + (Pi + Pj - 4u)p lPj 



3 \ u(u - l)(pi - Pj) 2 

+ Qw - (pj + pj){v + w) + 2p l p J )y{p i )y{p J ) \g u g u 
u{u-l){j) i -v)(p l -w){p J -v){pj-w){p l -p j ) 2 ) 1 3 

determine a three-field GT system. 

The systems in Examples [3] and 0] are related to the moduli spaces of curves 
of the respective genera g = 1 and g = 2. If we parameterize elliptic curves using 
the formula y 2 = x(x — l)(x — u) instead of the modular parameter, then we can 
rewrite the system in Example [3] in a form similar to the system in Example 21 

The functions / and h have poles at pi = pj in Examples HHU But there are 
GT systems holomorphic at pi = pj . 

Example 5 The formulas 

diPj = 0, diU m = g m (pi)diiii, did 3 ui=Q (37) 

define an n-field GT system for any n and for arbitrary functions g m (x). 
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3.2. GT-system extensions 

Definition 6 An extension of GT system ([25)1 using the field functions . . . , u n j rr 

is an auxiliary system of the form 

diu k = g k (pi,ui,...,u n+m )diUi, k = n + 1, . . . ,n + m, i = l,...,N, (38) 

such that Eqs. (|28[) and (|38[) are compatible. 



Example la (continuation of Example [l} . System (|3"3"|) admits the one- field 
extensions 

^ = v{v- 1) 
Pi - w 

a»« = (i + -W« ( 4 °) 

VPi Pi/ 

^^((^tf-^t)^' (41) 

Siti = f-^-r + - + (c 3 e" - c 4 e 2 > l + c 4 e 2 > 2 ) ^ ( 42 ) 

where Ci are arbitrary constants. 

Describing all extensions for a given GT system is a complicated problem. 
But extensions of type (f39|) always exist and can be constructed in a uniform 
way. Namely, it turns out that if we complete system (|28p with the equation 

diu n+ i = f(pi,u n +i,ui, . . . ,u n )diiii 

(the function / is the same as in (|28p). then the obtained (n + l)-field system 
remains compatible. We call this procedure the regular extension. For example, 
the regular extension in the case of Example [T] results in the additional equation 

a 2 v 2 + aiv + a 

OiV = OiU, 

Pi - v 

We note that if we take v as a new field variable, then the first formula in 
Example [1] becomes 

(a 2 p) + a 1 p ] +a )(p z - v) . . . , 

OiPj = ~, 1- ; T7 rOiU, i, J = 1, ■ • ■ ,N, l^J. 

(a 2 v z + a\V + a a )(pi - pj) 

We see that we can use this trick to obtain the canonical forms of the GT system 
in Example [5] from the corresponding GT systems in Example [TJ 

Clearly, we can apply the regular extension procedure repeatedly. For ex- 
ample, the n-fold regular extension of GT system © is 



— - — diW, didjW = — -7rd l wd j w, diu k = ® %W . 

Pi-Pj iPi—Pj) ' Pi-Uk 
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For this system, we can construct degenerations corresponding to merging some 
fields. In particular, merging all the fields Uj — > w, i = 1, . . . , n, results in the 
GT system 

dip, = — - — diWi, didjWi = — ^d l widjWi, 

I', I'. [Pi-PjY 

diWj = - ku k pj~ k ~ 2 j djWi . 

^ fc=i ' 

This extension at n — 1 exactly coincides with the GT system for dispersionless 
KP equation ([5]) (see Example [B] below) . 

We note that we cannot obtain the GT system in Example 2] from a system 
with fewer fields using a regular extension procedure. 



4 From GT systems to integrable models 

Constructing compatible families of systems of form (|M|) is our next step on the 
path from GT systems to integrable models (JXJ) . 

4.1. Families of (l+l)-dimensional hydrodynamic-type systems as- 
sociated with a GT system 

Definition 7 A family of semi-Hamiltonian systems of hydrodynamic type 

rt=F(p i ,u 1 ,...,u n )ri (43) 

parameterized by solutions piir 1 ,r N ) , i — 1,...,N, Ujir 1 , . . . , r N ), j = 
1, ...,n, of GT system (f28| is called a family of (l+l)-dimensional systems 
associated with system (f2"5)) . 

Two (l+l)-dimensional systems (jM]) are compatible if and only if the func- 
tions F and G satisfy functional equation (|25p . This equation has trivial solu- 
tions F and G = C\F + c%, where F is an arbitrary function and Ci are arbitrary 
constants. 

Not every GT system admits nontrivial solutions F and G of Eq. (|25|) de- 
pending essentially on the field functions ui. For example, such solutions exist 
only under a special choice of the functions g m in Example [5j 

We consider one-field GT systems (|3"TT) . In accordance with (|23)l . we must 
then have functions of two variables F(x,u) and G(x,u) such that 

f( x , y ,u)^(y,u) + ^(y,u) = f(x, y, u)jg (y, u) + u) 

F(x,u) - F(y,u) ~ G(x,u)-G(y,u) ' ( ' 

Proposition 1 Every one-field GT system admitting a nontrivial solution F(x, u), 
G{x,u) of Eq. (|4"4")) is equivalent to a system in either Example]^ or Example^ 



14 



The following results on constructing solutions of functional equation ([31)1 
for n- field systems with arbitrary n were obtained in [20 , [28 . 

Example lb (continuation of Example [1]). An n-fold application of the 
regular extension procedure to GT system (f55|) with a quadratic polynomial of 
general position results in the GT system 



dipj = — — - di-w, d l d j w = — ^-^-diwdjW, (45) 

Pi-Pj (Pi-Pj) 2 



i,j = l,...,N, i^j, 

d iUj = UjiUj ~ mW , i = l,...,N, j = l,...,n. (46) 

Pi - Uj 

We consider the compatible overdetermined system of partial differential equa- 
tions 

d 2 h s, dh Sk dh _ , , 

-, i,j = l,...,n, j^fc, (47) 



du.jduk Uj — Uk duk Uk — Uj du 



d 2 h _ ( ^ \ Sj Sj ^ u k {u k - 1) dh 



dujduj \ ) Uj(uj — 1) u j( u j — 1) j~f u k — Uj duk 

+ ( v Sk + Sj + Sn+1 + Sj + s »+ 2 ] dh (48) 

"■ U k Uj Uj ■ ~ 1 J DUj' 

k^j 

where 8\,..., s n +2 S C are arbitrary parameters. Solutions of this system be- 
long to the class of generalized hypergeometric functions |H] . We can easily see 
that the vector space H of all the solutions of this system has the dimension 
n + 1. For any h G H, we set 

n > n+2 v 

S(/i,p) = ^Mi(uj-l)(p-tti) (p-u n )h Uz + ( 1+^2 s i ) ■ ' ■ {p~u n )h, 

»=1 ^ i=l ' 

where i means that we omit the corresponding multiplier in the product. The 
function S is obviously a polynomial of degree n in p. 

Let /ii, foj and /13 be linearly independent functions in Then the func- 
tions 

_ S(fei,p) _ S(/t 2 ,p) v 

satisfy functional equation (|2"5")) . 

Moreover, we have solutions of Eq. (|25l) of form (|5"0")) for which S has a degree 
less than n. In this case, the polynomial S is to be defined in a more complicate 
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way. Namely, we fix k linearly independent functions hi, . . . , hk ETi- Then for 
each h E %, the polynomial S(h,p) is given by the formula 



S ( h ,p) = -T ^2 Ui(Ui - l)(p- Ul) ■ 



■i---(p-u n - k+ i)Ai(h), (51) 



where h G % and 



A = det 



/ 




h 2 




\ 




hi 


«„- fc+ 2 h 




hk,u n ^ 


fc + 2 




V hi <Un 






, / 






( k 


hi 










h Ui 




h2,ut 






det 










^■fc,ll„_fc + 2 








^2,ii„ 







We note that linear 
■ • + rffe/ifc with 



Clearly, S(h,p) is a polynomial of degree n — k in p 
transformations hi — > Cnhi + • • • + Cife/ifc and h — > h + d\hi 
constant coefficients and di do not change S(h,p). 

There are solutions of functional equation (|25|) for GT system (j45|) . 146[) 
different from those constructed above. For example, for n = 2, we have the 
solution 

0-mi)0-u 2 ) O-^OO-^) 



G 



p(u 1 -u 2 ) ' 0-l)(wi-u 2 )' 
Apparently, these solutions are related to singular rational algebraic curves. 

Example 3a (continuation of Example [3J . A regular rt-field extension of 
GT system (j35|) is given by the formulas 



d a u/3 = -—(p(p a - up) - p(p a ))d a T, /3 = 1,. 
2m 



(52) 



A nontrivial solution of functional equation (|25p is given by formulas (jST))) with 



(si + '" + s " ) wo) 



1G 



Here, 77 = SiUi+- ■ ■+s n u n +rT+r/o, where s\, . . . , s n , r, tjq are arbitrary constants 
and h(ui, . . . , u n , r) is a solution of the elliptic hypergeometric system 

K a up = sp{p{up - u a ) + p{u a + r)) - p(up) - p(r)))h Ua + 
+ s a (p(u a - up) + p(v,f3 +r))- p(u a ) - p{ri))h Up , 

K a u a =S a ^] (P(u a ) + P(V) - P( u a - Up) - p(up + l]))h Ufj + 



^ spp(u a - up) + (s a + l)p(u a + 77) + s a p(-rj) 



+ (sq - s a - l )p( u a) + ^irirj h Ua - s s a (p'(u a ) - p'(rj))h, 



We can show that the solution space of this system has the dimension n + 1. 

We construct elliptic analogues of solutions (l50l) . (fSTj) as follows. We fix k 
linearly independent solutions hi, . . . ,hk of the elliptic hypergeometric system. 
For any solution h, we then have 



' n — k 



A(»?)fl(p) 



(53) 



where 



A = dot 



/ h u 



A a (h) = dot 



( h 



A (h) = dot 



hi 



hi, u „ 



hk,u n - k + 1 
h k,u n I 

h k \ 

hk,u„ ) 



4.2. Integrable three-dimensional systems generated by a pair of 
compatible families. Given GT system ([28)) and a pair of compatible fami- 
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lies (pM]l (i.e., a solution F, G of functional equation ([25])). we can construct the 
corresponding system ([T]) by a simple calculation. 

Lemma We consider the linear space V of functions of the variable p generated 
by 3n functions 

{F(p,ui, . . . ,u n )gj(p,ui, . . . ,u n ), G(p,ui, . . . ,u n )gj(p,ui, . . . ,u n ), gj(p, m, . . . ,u n ); j = l,...,n}. 

Here, g\ = 1 by definition. The system of form Jl} with hydrodynamic reduc- 
tions (|24p then consists of n+k equations if and only if the dimension of V is 
2n — k. The coefficients of system (fTJ) are determined by the relations 

n 

^(aiji^Fip,^, . . . ,u n ) + b lJ (u)G(p,u 1 ,...,u n ) + c lj (u))g j (p,u 1 , . . . ,u n ) = 0, 

3=1 

(54) 

i = 1, . . . , n + k. 

Corollary // the functions F and G are determined by formulas of type (|50p . 
where the expression S(h,p) is linear in h and hi, and h% are arbitrary 
elements of a vector space, then the linear transformation hi — 2j=i kijhj 
corresponds to transformation (jTTJ) in the corresponding system (JTJ) . 

Example 6 We consider the extension diV — PidiU of GT system ([5]). The 
corresponding functional equation (|44|) has the solution F = p 2 + u, G = p. 
Moreover, we have g\ = 1 and g 2 — p. The space generated by the set of 
functions {p 2 + u, p 3 + up, p, p 2 , 1, p} has the dimension a — 4. The matrices 
A, B, and C composed of the unknown coefficients of system |T]) are determined 
up to simultaneous multiplication by an arbitrary matrix from the left. We can 
easily verify that using such a transformation, we can reduce any solution of 
the system of eight equations, which is equivalent to relations (|54p . to an = 1, 
a\2 = a 2 \ = a 2 2 = 0, bu = ^22 = 0, 612 = -1, &21 = 1, c\\ = -u, C12 = c 2 i = 0. 
and C22 = — 1. This solution corresponds to system (|5j). 
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An analogous calculation for ([45)1. (HSf . PD]) . and ()50l) results in the system 
for u\, . . ., u„: 

v-, - lKt, -m(ui-lKt, 

rr! Uj -Ui 

+ o(h q h r ^ Uj — h r h q , Uj )Ujf t + 

Uj(uj - l)t*x,t g -Ui(Ui - l)Uj,t g 



+ a{h r h s ^ Uj - h s h r ^ Uj )u 0ttq + 

U 3 (Uj - l)Uij r - Ui(ui - l)Ujj r 



+ a(h s h qiUj - h q h s . u .)u %tr = 0, (55) 



where j = 1, . . . , n and a = 1 + s\ + ■ ■ ■ + s n+ 2- Here, the subscript after 
the comma indicates the partial derivative with respect to the corresponding 
argument, q — 0, r = 1, s = 2 and to = x, t\ = t, fa = y in expressions like h qiUj 
and Ui, tq ■ 

Remark Let ho, hi, hi, . . . , h n be a basis in the space of solutions of hypergeo- 
metric system (|47|) , (|48p . We associate the proper time ti with each hi . For each 
triple of pairwise distinct < q,r, s < n, formula (1551) defines the corresponding 
three-dimensional system with the times t q , t r , and t s . According to I20f , all 
these systems are mutually compatible. 

The explicit form of systems related to the case k > can be found in [2D] . 
Systems generated by the elliptic GT system in Example[3]were presented in [28] . 

5 Integrable weakly nonlinear (2+l)-dimensional 
systems 

The functions / and h have a pole on the diagonal Pi = pj in all GT-system 
examples (|2"5)) above. But we also have GT systems holomorphic at p, t = pj. In 
particular, such are the systems in Example [5] 

Definition 8 Integrable system ([I} is said to be weakly nonlinear if the corre- 
sponding GT system is holomorphic on the diagonal pi = pj. 

The GT-system holomorphicity conditions are expressed as follows in terms 
of the solution F, G of functional equation ([25]) . Writing (|25l) in the form 



(d iPj F'( Pj ) + F( Pj ) Ua d lUc }j (G( Pt ) - Gfa)) = U Pj G'( Pj ) + G( Pj ) Ua d lUc }j (F( Pt ) - F( Pj )), 
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where the prime denotes the derivative with respect to pj , and expressing diPj 
from it, we find that this expression is holomorphic on the diagonal if and only 
if 

Y,{F(Pi)u a G'( Pl ) - Gipi^F'ip^diUa = 0. (56) 

Proposition 2 Let system Q]) be weakly nonlinear, and let k — (i.e., the 
number of unknowns is equal to the number of equations in the system). Then 
the matrix 

A = (a n A + a 12 B + ai 3 C) _1 (a 2 iA + a 22 B + a 23 C) 

is linearly degenerate for all constant aij . 

We recall that a matrix P(u\, . . . , u n ) is said to be linearly degenerate if the 
derivative of any eigenvalue along the corresponding eigenvector vanishes. 

Proof Using substitution (j29[) . we reduce G(p) to the form G = p. Then (|56p 
becomes 

^F( Pi ) Ua d iUa = 0. (57) 

a 

An arbitrary change of variables (jlip brings system |T]) to the form u t = Au x + 
Bu y , where 

A = (a n A + a 12 B + ai 3 C) _1 (a 21 A + a 22 B + a 23 C), 

B = (a n A + a 12 B + a 13 C)- 1 {a 31 A + a 32 B + a 33 C). 

For fi l — pi and X 1 = F(pi), formula (|21[) for this system yields 

(A + piB - F{ Pl ) Id)(ftwi, . . . , d t u n ) T = 0. 

For this system, formula (|57[) implies that the derivative of the eigenvalue F(pi) 
along the eigenvector (diU\, . . . , diU n ) vanishes. The matrix A+piB is therefore 
linearly degenerate at any value of the parameter pi . 

The obtained result states that if the number of equations in system ([1} is 
equal to the number of unknown functions, then every two-dimensional system 
describing solutions u = u.(c\x + c 2 y + c 3 t, c^x + c 5 y + c^t) of the traveling-wave 
type for weakly nonlinear system (jTJ) is weakly nonlinear in the sense in |40j . 

Example 7 We consider the two-component system (see J13f) 

{wtj + (o bj {w^j + i^r sj (wyj °' ^ 
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where 

P(w) P(v) 

a = to, b = v, r — , q = , 

w — v v — w 

s= h-P [v), p= \--P{w). 

w — v 6 v — w 6 

Here, P is an arbitrary third- degree polynomial. The corresponding GT system 
P(w) , , ( 1 , P'(v) 



(w — v)P(v) J \w — v P{v) 
1 P'(w)\ P(v) 



v — w ' P(w) J J (v — to)P(to) ' 

<9i« = p;d;to, 

/ P(w) 1 P'W ^ a 

djdjW = - . . . piPj H h , , )OiwdjW 

J \ (t? - to)P(t;) u-to P(to) / 

is polynomial in Pi andpj, and the initial (2+1) -dimensional system is therefore 
weakly nonlinear. We can verify that this GT system is equivalent to the GT 
system 

diPj = 0, diU 2 =PidiU\, didjiii = 0, 
which belongs to the class in Example^ 

As shown in [T3] , system (JSJ has a linear pseudopotential representation with 
the spectral parameter. This turns out to be a general property of systems (|3T[) 
in Example [5j 

Proposition 3 Let F(p, u±, . . . , u n ) and G(p, u±, . . . , u n ) be a solution of func- 
tional equation (|25[) for a GT system of form (|37|) . Then the corresponding 
(2+1) -dimensional system admits the pseudopotential representation 

tpt = F(X,ui,...,u n )ip x , ip y = G(X,ui,...,u n )ipx, (59) 
where X is the spectral parameter. 

Proof From (|25|) with the vanishing of diPj taken into account, we obtain 
J2 a F (P])u a diU a _ J2 a G (Pj)^A u a 



F{ Pi )-F{ Pj ) G( Pt )~G( Pj ) 



(60) 



On the other hand, the compatibility condition for system (|59p is F y + FG X 
Gt + GF X , whence we obtain 

( F ( X )^ u ay + F(X)G(X) Uot u ax ) = Y, ( G ( A )- Q u at + G(X)F{X) Ua u ax ) . 
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Substituting the expressions 

u ay = ^ d. t u a r l y = ^ diU a G(pi)r l x , 



u at = ^2 d l u a F(p i )r l x 

i i i 

in this equality and replacing A with pj , we obtain (|60p . 

5.1 The general position case 

Using our observation on the equivalence between GT system (|58[) and a system 
of form (1371) with rational functions g m , we now generalize Example [7] to the 
case of arbitrary n and k. 

We consider an (n+l)-ficld GT system 



dipj = 0, diu m = — — — -diW, didjw = (). 



Pi - A, 



(61) 



where Ui,...,u n and w are the field functions and Ao, Ai, . . . , X n are the pa- 
rameters. We let H n denote the linear space of functions in the variables 
ui,...,u n generated by the elements 1, e Ul , . . . , e Un , For any function g — 



a + aie 



a n e Uri S H n , we set 



S n {g,p) 



P - A p - A 



E 



For any k € N such that < k < n — 1, we fix the functions hi, . . . , hk & H n , 
where hi = bio + bi ie Vl + • • • +». n e u ™ , and set 



S n ,k{g,p) = dct 



(S n {g,p) S n (hi,p) S n (h 2 ,p) 
g hi h 2 

a n-k+2 &l,n-fc+2 &2,n-fc+2 



By definition, S nfi (g,p) = S n (g,p). 



Sn(h k ,p)\ 
hk 

°k,n-k+2 
bk,n j 



(62) 



Theorem Let gi, g 2 , and g% be linearly independent functions in T-L n . Then 
for any < k < n — 1, the functions 



F 



S n ,k{gi,Pi) 
S n ,k(93,Pi) ' 

satisfy functional equation (|25l) . 



G 



S n ,k{92,Pi) 

Sn,k(g3,Pi) 



(63) 



We now write the corresponding (2+l)-dimensional systems explicitly. Ac- 
cording to the lemma in Sec. 4.2, the equation 



^2(AiUi ttl + BiU i<t:i + CiUi >x ) = 
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belongs to the system if and only if the expression 
" A l - A 

V" — — 7-{AiS n! k(gi,p) + BiS n! k(g2,p) + CiS n . k (g 3 ,p)) 

vanishes identically. We first consider the case k — 0. Let gi — <ij,o + ai ie Ul + 
■ • ' + a i,n eUn i i = 1) 2, 3. Then the corresponding (2+l)-dimensional system is 

n 

E, \ m U i,ti ~ u j,ti , , \ M i,ti 
(02,103,3 - a 2 ,ja 3 ^)e 3 — h (a 2) ia 3i0 - a 3ii a 2 ,ojT r- + 
Ai — Aj Ai — An 

3=1, 

n 

Uj U i,t2 ~ U 3, t 2 I l \ U i,t2 | 



- y{a.3.iaij - as.jai^je 3 — \- [aa^aifi — ai,ia3,o) 

rr — A 7 Ai — An 



+ V(«i,i«2,i - aija 2 ,i)e 3 — h (ai,,a2,o - a 2 ,iai,o)-r r- = 0, 

z — ' A; — Aj Ai — An 



3=1, 



(64) 



where i = 1, . . . , n. In the simplest case, n = 2, and we can use linear transfor- 
mations of the functions gi (see the corollary in Sec. 4.2) to reduce the matrix 
ciij to the diagonal form diag((Ai — Ao)" 1 , (A — A2) -1 , (A 2 — Ai)" 1 ). System (j64|) 
then takes the simple form 

u y = e v (v x — u x ), v t — e u (u x - v x ). 

The system in Example[7]reduces to this system upon some changes of form (| 1 1[) . ()10[) . 
Introducing the potential 

ry U ry V 

Zt = e, Z y = -e , 
we can write this system as the single equation 

Zty Z X y Z x i 



ZfZy Zy Zt 

If k > 0, then the corresponding (2+l)-dimensional system contains two 
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groups of equations: 

n-k+l 

£ (A 4 ( 52 )A,( 53 ) - A 3 -( g2 )A i ( fl3 ))e^ U< -; i ~" 3 - tl + 
+ (A 4 ( 32 )A (.g 3 ) - A (g 2 )A,( 33 )) 



H,ti 



n-k+1 



+ (A J ( ff3 )A J ( ffl )-A J ( 53 )A,( 5l )) e " 



. Uj tt2 - u jtt2 



3 = 1; 



n-k+l 



Aj — A, 



+ (A i ( ff3 )A ( 5 i) - Ao(ff 3 )Ai(fli)) 



A,-A 



£ (A,( 5l )A J ( 32 ) - A J -( 5l )A i ( 92 ))e«^ 



3,x 



3 = 1, 



A, - A< 



+ 



+ 



+ 



+ 



+ ^( 5l )Ao(.9 2 ) - A (g 1 )A i (g 2 ) S j = 0, 



where i = 1, . . . , n — fc, and 

n— fc+1 n— fc+1 

j=i i=i 



ra-fc+1 



n-k+l 



Aj — Aj 



A, — Aq 



Uj,t r ~ Uj,t r 

Xi - Xi 



+ A 



where i = n — k + 2, . . . , n. Here, r, s — 1, 2, 3, t 3 = x, and 



Aj(ff) = det 



/ 5 ^1 h 2 

a n-k+2 t>\,n-k+2 ^2,n-fe+2 



.9 = a + aie'" 1 H h a n e Un , 

hi = &i,o + + • • • + bi,ne u " 



/ifc \ 
bk,j 

^fc,n-fc+2 



h k = b kfi + b kil e Ul H h &fc,„e u " 



The first group contains n— k equations, and the second group contains 3k equa- 
tions. There are exactly k linear combinations of equations in the second group 
that contain no derivatives of the Ui, i = n — k + 1, . . . , n. These linear combina- 
tions can be linearly expressed in terms of the equations in the first group. The 
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(2+l)-dimensional system therefore contains (n — k) + 3k — k = n + k linearly 
independent equations. 

The simplest nontrivial case is where n = 3 and k = 1 . In Appendix ^ we 
show that the corresponding (2+l)-dimensional system is equivalent to Eq. ([7]). 



5.2 Degeneracies 

Gibbons-Tsarev system (fuTl) and formulas (p3")l correspond to the case of pair- 
wise distinct parameters Aq , . . . , A n . Here, we consider some reductions of the 
construction in Sec. 15.11 that correspond to merging these parameters. 

We define the polynomials Pi(ui, U2, ■ . . ) as the coefficients of the Taylor 
expansion 

e eu 1+ e 2 u 2 + - = l + p i£ + p 2£ 2 + _ 

In particular, 

P 1 (u 1 ,u 2 , ...)=ui, P2{ui,u 2 , ■■■) = Pz{ui,u 2 , ■■■) = u z +uiu 2 +^u\. 

We let <3fc(e, ui,u 2 ,. . .) denote the partial sums 

1 + Ei=i Pi**- We set P Q = 

Qo = 1 by definition. 

Let Ao , . . . , A; be pairwise distinct roots with the respective multiplicities 
n + 1, rii, ■ • ■ i n i- Clearly, n + • • • + n\ = n. The corresponding (n+l)-ficld GT 
system with the field functions «o,i) • • • > u o,n a , Ui l, . . . , u% >ni , . . . , Ui t %, . . . , Ui iHl , w 
is ^ 

d t pj = 0, diU ,m = t 7-^—diW, 

[Pi - A ) m 

(65) 



diU s m — t — h t — - t diw, didjw = 0. 

In the case of GT system (p5|) , the solutions F, G of functional equation (|25|) 
are given by the formulas in the theorem above where the functions S n ^ are 
defined as follows. We let H denote the vector space generated by the functions 
P<(iio,i,uo,2,...)) i = 0,...,no, and e u "' 1 Pi- 1 (u a<1 ,u Sj2 , . . . ), s = i = 

1, . . . , n s . For any element 

no / n s 

g = ^2ao,iPi(u ,i,u ,2, ■■■) + ^^a s , i e Ms ' 1 P i _i(u^i,u Si2 , ...) 

i=0 s=l i=l 

in this space, we set 

a I \ ST" a 0,iQi{P~ Ao,tiO,l; u O,2,...) , a^e" 3 ' 1 Qj_i (p - A s , U al , U S)2 , . . . ) 
S n {g,p) = > -, r-^-TT h > > 7 r-r^ . 

U (p- x oY +1 (p-A.)* 

By definition, S nfi (g,p) = S n (g,p). 
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In the case k > 0, we fix the elements 



no 



l n s 



hj = ^&0,i-Pi(U0,l,U0,2, • ■ ■)+^2z2 b s,i eUs ' 1 Pi-l(u s ,l,U s ,2, ■ ■ •)> j = 1, • • • , k, 



i=0 

in Then 



(S n {g,p) S n (hx,p) S n (h 2 ,p) 
g hi h 2 



Sn,k(9,P) = dct 



V„_fc +2 



S n (hk,p)\ 

hk 



(66) 



where (vi,...,v n ) = (u ,i, . . . , Wo,n OJ "i,i> • • • > • • • : • • • > We can 

then easily reconstruct the explicit form of the corresponding (2+l)-dimensional 
systems using the general construction in the lemma (see Sec. 4.2). We omit 
this answer for brevity. 



6 Discussion: Classifying integrable (2+l)-dimensional 
systems of hydrodynamic type 

The path from GT systems to integrable (2+l)-dimensional systems of hydro- 
dynamic type outlined in Sees. [3] and 2] can be regarded as the following project 
of classifying integrable (2+l)-dimensional systems. 

Step 1 Classify all GT systems (|28[) (up to equivalence relations (|29p and (|30[) ). 

Step 2 Classify all the nontrivial solutions F,G of functional equation (|25[) for 
each GT system. 

Step 3 Use the lemma construction to construct the (2+1) -dimensional systems 
corresponding to each GT system and each solution of Eq. (|25|) . 



In the following remarks, we share our experience acquired during several 
years of working with integrable systems ([1]). The major part of these notes are 
observations and conjectures. We currently have few rigorously proved state- 
ments. Most of unproved statements are difficult. In particular, the problem 
of classifying n-field GT systems with arbitrary n is very difficult. But such a 
classification is apparently possible for lower n. We intend to study this problem 
in the future. 

Remarks concerning Step Q] 

1. The different general-position GT systems are few and are universal: sev- 
eral families of (2+l)-dimensional systems, each family depending on essential 
parameters, correspond to each of the GT systems. The known GT systems 
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correspond to algebraic curves. Apparently, for any genus g = 0, 1, . . . , we have 
a unique GT system with one field for g = 0, 1 and 3g — 3 fields for g > 1 (see 
Examples HHU such that every integrable general-position (2+l)-dimensional 
system corresponds to either this GT system or its regular extension. From 
the algebraic-geometric standpoint, a regular extension corresponds to adding 
a marked point on the curve. 

2. The definition of GT system using formulas (|2"5)l admits a coordinate- free 
reformulation. Let M be a fiber bundle with a one-dimensional fiber E and an 
n-dimensional base F . Every pi is then a coordinate on E, and u\, . . . ,u n are 
the coordinates on F. We can thus define the notion of the GT structure on M . 
We then apparently have a canonical GT structure on the natural fiber bundle 
over the moduli space M g of genus- <? curves. A fiber E of this fiber bundle is 
the curve corresponding to a point in M g . 

3. Integrable systems |T]) are defined up to arbitrary pointwise transfor- 
mations u — > $(u). The equivalence problem for systems ([T]) is nontrivial; in 
particular, finding a coordinate system in which a (2+l)-dimensional system 
has the simplest form is often difficult. One possible approach to this prob- 
lem is based on the observation that the simplest coordinates for a GT system, 
as a rule, are simultaneously the simplest ones for the corresponding (2+1)- 
dimensional system. 

4. "Proper" GT system coordinates admit an algebraic-geometric interpre- 
tation. For example, we can interpret each of the variables pt and u in Example[2] 
as a coordinate on CP 1 and each of the variables pi in Example[3]as a coordinate 
on an elliptic curve with r being the modular parameter of this curve (i.e., a 
coordinate on the moduli space of genus-one curves). The field functions u, v, 
and w in Example |4] are coordinates on the moduli space of genus-two curves. 
We expect that in the case g > 2, we have a GT system whose 3g— 3 field 
functions are coordinates on the moduli space M g of genus-g curves. Finding 
good formulas here is hindered by the absence of a good choice of coordinates 
on M g . We can obtain an algebraic-geometric description of these GT systems 
simultaneously with the description of hydrodynamic reductions of the (2+1)- 
dimcnsional systems from Krichever's paper [6]. As soon as we find GT systems 
related to curves of genus g > 2, describing solutions of the corresponding func- 
tional equation (|25p would apparently result in constructing (2+l)-dimensional 
systems depending on essential parameters. We must then obtain the first ex- 
amples of generalized hypergeometric functions related to higher-genus curves. 
This program was realized in [20], [28] for g = 0,1. We note that the (2+1)- 
dimensional systems in [6] do not contain essential parameters and correspond 
to the simplest integer values of parameters in the hypergeometric series. 

5. The problem of describing all the GT systems with a given number of 
fields includes the problem of describing degenerations of a general position 
system. Degenerations appearing when merging the field functions Ui (merging 
the marked points in the algebraic-geometric language) can be described using 
Young diagrams. In Sec. 3.2, we considered the degeneration of the GT system 
in Example [2] under which all the points merge into a single point. Gibbons- 
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Tsarev system (|oTj) corresponds to degeneration of the type m — > Aj + eui, 
e -> 0. 

Remarks concerning Step [2] 

6. Algebraic curves in the classification of integrable (2+l)-dimensional sys- 
tems appear in Step[5J In all the known examples, solutions F, G of functional 
equation (|25p are algebraically dependent as functions of p. The polynomial 
relation between them is called the dispersion curve (see (|2"2")) ). If we are inter- 
ested in solutions F, G parameterizing a curve of a given genus for a given GT 
system, then we assume that the dependence of F, G on p is known. For exam- 
ple, in the case of a rational curve, the functions F, G are rational functions of 
p. Different choices of the degrees of F, G result in different families of (2+1)- 
dimensional systems. The functional equation becomes a system of nonlinear 
differential equations for the ansatz coefficients after the corresponding ansatz 
is substituted. 

7. If we additionally assume that system (Q]) admits a pseudopotcntial repre- 
sentation, then we can find solutions F, G using functional equation (|76p , which 
is simpler than (|25l) (see Appendix IB)) . 

8. The class of integrable systems ([l]) is closed under changes of vari- 
ables (fTTj) . The existence of this group leads to formulas ([50| for solutions 
F,G of functional equation (|2"5j) . In these formulas, the expression for S(h,p) 
is linear in h, and hi, h2, and /13 are arbitrary elements of a vector space Ti of 
dimension M. In the general position case, 1-L is the solution space for an overde- 
tcrmincd system of linear partial differential equations of hypergeometric type. 
We can interpret formulas (|50")l as a linearizing substitution for the nonlinear 
system mentioned in item 5. The space H may consist of elementary functions 
in the case of reduction. For example, the space H consists of quasipolynomials 
in Sec. [5j 

9. In the case of curves of genus g = 0,1, the dimension of H is n + 1, 
where n is the number of marked points. The set of times tj, i = 1, . . . , n + 1, 
in H corresponds to the basis hi, i = 1, . . . ,n + 1. For any triple tj, tj, tk of 
times, formula ([50)1 results in a three-dimensional system of form ([1]) with the 
independent variables ti, tj, and tk- All the systems corresponding to different 
triples of times are mutually compatible. 

A Conservation laws of hydrodynamic type and 
integrable quasilinear second-order equations 

A relation of the form 

a(u) t + b(u) y + c(u) x = 0, 

where a, b, and c are some scalar functions, that is satisfied by any solution u of 
system (|55[) is called a conservation law of hydrodynamic type for system (JTJ) . 
All known integrable systems ([1} have a rich collection of conservation laws. 
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More precisely, we consider systems ([T]) with k = 0. It was proved in [2UJ 
that systems ([55l) corresponding to a rational nonsingular curve (Example [lj 
have n+1 linearly independent conservation laws; all partial derivatives of the 
functions a, b, and c with respect to the it, were found explicitly there. The 
dimension of the space of conservation laws is apparently n in the case of an 
elliptic curve (Example [3]). This was verified for degenerations of the elliptic 
systems presented in [7|. 

In the case of exponential systems (|64[) with k = 0, we can find conservation 
laws explicitly. 



Proposition 4 System (|64p with k = admits the n+l conservation laws 



awe -1 *' a 2 .ie~ Mi 



where i = 1, 



, n, and 



O2,0 ^3,0 



a 3,0 a l,0 



Ol,0 a 2,0 



Here, 



^(g,X i )=[S n (g,p) 



P - A* 



5; eg ( ff ,A )= 5„( fl ,p)- 



rt 



p - A 



p=A 



In the case fc > 0, all known integrable systems ((TJ) have fc triples of "short" 
conservation laws of the form 



a t — c a , 



We can introduce the potential Z x — a, Z y — b, and Z t = c for each such 
triple. For small k (for n > 2k in the rational case and for n > 2k — 1 in the 
elliptic case), introducing potentials Z%, . . . , Zk, we make system ([T]) a properly 
defined system consisting of several second-order equations and several first- 
order equations. In particular, for n = 3fc, system ((T]) is equivalent to a system 
of the form 



AiZtt + A 2 Z yt + A 3 Z xt + A A Z y 



A§Z X y 



A&Z XX — 0, 



(67) 



where Z = (Z\, . . . , Zk) T and Ai are some kxk matrices depending on Z x , Z y , 
and Z t . 

Short conservation laws for the systems in Sec. 4.1 with k > are given by 
the formulas 



( A(g r ,h 1 ,...,i,...,h k ) 



V A(/n 



,h k ) 



A(g s ,hi,...,i,...,h k ) 
A(hu...,h k ) 



(68) 
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where i = 1, . . . , k, r, s = 1, 2, 3, t\ = t, ti = y, ta = x, and 



A(/ 1 ,...,/ fe )=det 



/ h 

V fl,u. 



fk \ 



h,u„ 



h 



fk,i 



fk 



J 



In addition to these conservation laws, such systems apparently also have n+1 
linearly independent conservation laws of hydrodynamic type. 

In the case where n = 3 and k = 1, the system in Sec. 4.1 depends on four 
arbitrary linearly independent constant vectors (see ([511)). After the potential is 
introduced using the short conservation laws, the system becomes an equation 
of form ([2]) whose coefficients Ai are polynomials of a degree not exceeding two 
in the variables Z x , Z y , and Z t . Using the group GL(4) acting on the set of 
equations @, we can reduce the matrix composed from the above vectors to 
the unit matrix by linear transformations of the variables Z ', x, t, and y. The 
most convenient way to do this is to make the function h equal to unity. Then, 
in the general position case where no = and n\ = 112 = TI3 = 1 (see (J65|), 
the equation coincides with ([7]). The degeneration of the form no = 3 and 
n\ = 712 = n3 = results in the equation [41) 

Z X t ~ Zyy + Z X Z X y ~ Z y Z XX = 0, 

the degeneration no = 2, n± = 1, and = nj, — gives 

Z y t + ZfZ xx — Z x Z x t — 0, 
the degeneration no = 1, n<i = 2, and ri2 = = gives 

ZyZ X t — ZfZ X y + Zyy = 0, 

and the degeneration no — ni = n-i — 1 and n3 = gives 

ZyZ x t — ZtZ X y + = 0. 



B The GT systems and the pseudopotential rep- 
resentation 

In Sec.[4j we described several classes of solutions of functional equation ([251 for 
the rational and elliptic GT systems. These solutions were found in [20], [28] un- 
der the additional assumption that in addition to hydrodynamic reductions, the 
corresponding (2+l)-dimensional systems |T]) have the pseudopotential repre- 
sentation or, equivalently, the dispersionless Lax representation. As we demon- 
strate below, this assumption results in a functional equation that is simpler 
than ([25]) and was in fact solved. 
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The dispersionless Lax representation is the relation 

L t = {L,A}, (69) 

where {L,A} = A\L X — A X L\. Here, A = A(X, u\{x, t), . . . , u n (x, t)) and the 
unknown function L = L(X, U\(x, t), . . . , u n (x, t)) depend on the spectral param- 
eter A. The transformation L(x, t, A) — > A(x, t, L) takes (|6"9"]) to the conservative 
form 

A t = A(X, ui, . . . ,u n )x, (70) 

where L now plays the role of a hidden parameter. We can rewrite the last 
equation as 

$ t = A($ x ,Ui,...,Un), (71) 

where A = $> x . 

We assume that system (JTJ) has both hydrodynamic reductions (l24l) and 
pseudopotential representation Calculating L t and L K on the iV-phase 

solutions and using ((33]), we write as 

ui, . . . ,u n )—Lx ^ Ai fc (A, mi, ... , u n )diUk \ i 

i i ^ fe=l ' 

Equating the coefficients of r x , we obtain 

n 

F(pi,ux, ■ ■ .,u n )diL = Ax(X,Ui, . . .,u n )diL - L\ A Uk (A, u±, . . .,u n )diU k . 

fe=i 

The function L therefore satisfies the system of equations 

d i L = t(p i ,X,u 1 ,...,u n )Lxd i u 1 , i = l,...,N, (72) 

which is compatible by virtue of the corresponding GT system (|2"g)l . System (j?2"j) 
is called the Lowner equation. 

Using transformation (f2U|). we can arbitrarily gauge the function F. For 
example, we can set F = Pi. It is convenient to choose 

Ffc, ui, ... , u„) = A Pi (pi, u x , . . . , u n ), 

and system (f72j) then becomes 

ft£= n=i^^^--;'»")^ r A , £ = l,...,iV. (73) 
Aa(A, ui, . . . ,u„) - A Pi {pi,ui, . . .,u n ) 

The transformation L(x,t, A) — > X(x,t,L) brings (|73p to the form 

*A= n = i^(A,"i, ; -,«n)^ i = 1 ,... tJV . (74) 

-<4-Pi lPi> «1> • • • ) Un) - 4\(A, Ul, . . . , U n ) 
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Relation (fT3"|) implies that L\ = at the points A = pi, i — 1,...,N. 
Substituting A = pi in (|69l) , we obtain 



(L(pu r\ . . . , r N )) t =A Pi (p i ,u 1 ,... i u n )(L(pi, r\ . . . , r^))^. 

We see that the functions L(pi,r 1 , . . . , r ) and satisfy system (Pf3")l with the 
same functions Fj. Hence, L(pi,r 1 , . . . ,r N ) — Hi(r l ) for some functions /Xj. 
After changing the Riemann invariants ^i{r l ) — > r l , we can set 

r l =L(p u r\...,r N ). 

After the transformation L(x, t,A) — > X(x,t,L), this relation becomes 

^ = A(a;,t, r 4 ). 

Substituting L = pj in ([71|) . we obtain the relation 

dv - ELl^fe^b-^n)^ .. . . ._ j „ 

^ p ,(k,ui, ...,u n )- Apjlpj,^, ...,u n ) 

(75) 

which expresses the function f(pi,Pj, iti, . . . , it n ) in the corresponding GT sys- 
tem (|28p in terms of the potential A and the functions 52, • ■ • , <7n- This formula 
for n = 1 is contained in |42j . 

We recall that the GT-system gauge is fixed by the condition F = A Pi in 
formula (|75[) . This gauge is generally inconvenient. After an arbitrary transfor- 
mation (I2U1) , formula (f75|) becomes 

g _ Efc=l( $ Pi(Pi> M l: ■ • ■ ,Un)4>u k (Pj,Ui, ...,U n ) - ^piiphUt, . . .,U n )$ Ut (j>j,«i, . . . ,U„))cW 
1 j <P Pi {Pi,Ui, . . .,U n )$ pj (pj,Ui, ...,U n )~ <t> Pj (Pj,Ui, . . .,U n )$ Pi (pi,U-L, ...,«„) 

(76) 

where i ^ j, i, j = 1, . . . , N, and the potential is defined parametrically 

A = 0(p,«i, • . . ,«„), A = $(p,«i, . . . ,w„). 

Using (|28j) and omitting the arguments iti,...,u„ for brevity, we rewrite 
formula (|76[) in the form 

,, \ = YTk=l < ^P^(Pi) ( t>u k {P]) - 4> P APi)®u k {Pj))gk{Pi) (77) 

where gi = l. For a fixed GT system (|28p. this relation is a functional equation 
for the functions and $ that is much simpler than functional equation (|25l) 
for the functions F and G. In particular, fixing pj, we can find the dependence 
of the ratio 4> Pi /$ Pi on the variable pi from (|77p. Every solution of Eq. (J77J 
determines an integrable potential corresponding to system (|28[) . 

Example 2a (continuation of Example [2J . We find the one-field integrable 
potentials corresponding to one-field GT system with P(x) = x(x — 1). 
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Fixing pj in ([77)1 , we find that 4> Pi / <J> Pi is the ratio of polynomials of the first 
degree in pi. We set 

$ Pi = (Ai(u)pi + A (u))Z(pi,u), <f> P( = (B x (u)pi + B (u))Z(pi,u). 

Substituting these expressions in (|77[) and equating the coefficients of like powers 
of pj , we find that ([77)1 is equivalent to the relations 

$ u = -(A 1 (u)u + A Q (u)) Pt{ ^ ~ X) Z( Pi ,u), 

u(u — 1) 

K = -(B 1 (u)u + B (u) f i ^ i ~ Z(p h u). 

u(u — 1) 

We find both partial derivatives of the function Z(j>i, u) from the compatibility 
condition for these two systems (i.e., from the condition that the mixed deriva- 
tives of 4> and <£> are equal) . We then find the condition for the equality of the 
mixed derivatives of Z . Equating the coefficients of like powers of p; in this 
relation, we obtain a system of nonlinear ODEs for the functions Ai and Bi, 
i = 0,1. We can express a solution of this system in terms of two arbitrary 
solutions yi and y 2 of the standard hypergeomctric equation 

u(u - l)y{u)" + [{a + (3 + l)u - j}y(u)' + a(3y{u) = 

as 

Ax = -ay!, A = u^-V^y'x+auyx, B x = -ay 2 , B = u(u-l)y' 2 +auy 2 . 
We reconstruct the function Z from its partial derivatives: 
Z(p, u) = (p - l) a+ ^p-^-^~\p - uf' 1 . 

We finally obtain 

*(p,«) = f P [u{u - l)y[(u) + a(u - Ol/i («)](£ - l)^-^ 7 "^ - u)" -1 d£. 
Jo 

(78) 

To obtain the function cf>, we merely substitute y 2 for y\ in this formula. 

The potential found above corresponds to the general position case in [27] . 
where all the integrable one-field potentials were described. The generalization 
of formula {7TJ to the case of n- field GT system (gSJ, (gBJ is 

p„(p,/o= r^o^-ui)"* 1 " 1 ---^-^)"""" 1 ^ 1 - 1 ^-!)-""^- 1 ^. 



where /i(ui, . . . , u n ) is an arbitrary solution of system (jT7| . P5|) and the poly- 
nomial S*„ is given by formula (|49|) . The elliptic analogue of this potential is 
given by 
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where we use the notation in Example [3J Integrable potentials with the de- 
fect index k > corresponding to rational and elliptic curves were presented 
in [3D], HE]- 

Clearly, for any solution 0, $ of Eq. (|77|) . the pair ki<f) + fc 2 $, fc 3 + &4$ is 
a solution for any constant ki. In the above example, such a transformation 
corresponds to another choice of the hypergeometric equation solutions 2/1,2/2- 
The function <fi therefore always lies in a two-dimensional space V. If this 
space has the dimension d > 3, then choosing linearly independent elements 
(j>i,<p2,<j>3 G V, we obtain the functions 

p = (Mm q = (h)vi 

{<f>3)pi ' (0s)p« 

satisfying functional equation (|25|) . This is how solutions of Eq. ([25} were found 
in [20], [28]. Equation ([25} for GT system ((6lJ) was solved directly. 



C Multidimensional integrable systems of hy- 
drodynamic type 

Multidimensional generalizations of systems ([!} are systems of the form 

£^(u)J^=0, d>3, (79) 

i— 1 

where Ai are nxn matrices and u = (ui,...,u n ). We do not here discuss 
the definition of the integrability of such systems based on the hydrodynamic 
reduction method (see [13). But if we know what integrable system (TTJ) is, we 
can naturally require that every reduction 

, d d d \ 

u = u(^fci jl x i , ^2k 2 .iXi, k3,iXj J (80) 

^i=l i=l i=l ' 

of integrable system (jT9"]l result in a three-dimensional integrable system (TTJ) 
for any constants kij. We note that the coefficients of this system of form (fTJ) 
depend on the parameters fcjj linearly. It was shown in |43j that in the case 
n = 2, systems (|79p have this property, being reductions of the system 

u Xl + v X2 + uv X3 - vu X3 = 0, u Xi + v X5 + uv X6 - vu X6 = 0. (81) 

System (|81[) admits a maximum collection of hydrodynamic reductions and 
can be written as the commutation condition for a pair of vector fields. This 
system was first derived in [15] as follows. We have four parameters (namely, 
the coefficients po, pi, P2, and p$ of the polynomial P) in system ([55} entering 
the system coefficients linearly. If we write ([58} as 

3 

u t + A(u)u x + y^pjBj^Uy = 

i=0 
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and pass to the six-dimensional system 

3 



u t + A(u)u x + Bi(u)u yi = 0, 



then this last system is equivalent to system (f8Tj) . 

Analyzing system ([5"8"]) and other known examples of multidimensional inte- 
grable dispersionless systems, we proposed the conjecture that for any system ([T} 
obtained by reduction (|80[) from an integrablc multidimensional system, the cor- 
responding GT system is equivalent to system (|3T)1 [H]. 

In Sec. [SJ we constructed a broad class of systems (|T|) whose corresponding 
GT systems have form (13"T1) . Systems (|p^|) are the simplest ones. Below, we at- 
tempt to construct multidimensional integrable systems corresponding to these 
systems. 

System (164[) is written in the coordinates Ui in which the corresponding GT 
system acquires the simplest form. According to Proposition [3J system ([Ml) 
admits a pseudopotential representation depending rationally on the spectral 
parameter. Taking zeros of the pseudopotential representation denominator 
as the new field functions Vi, i — 1, . . . , n, we reduce the system to the form 
analogous to ([58l) 

(vj)t- — i)=e+ ( ^ 7ikP(vk)+Si ) My+Y f3ijP(vi)(vj)y =0, i= 1, ...,n, 

where /3y, 7^, and (5^ are expressed in terms of the two symmetric polynomials 

/( v ) = /o + /i("H 1- V n ) H 1- fnVi ■■ - v n , 

ff( v ) = .9o +9i(vi H h u n ) H 1- 9nVi - --Vn 

(which are linear in each of the variables) and the polynomial P{x) = x n+1 + 
p n x n H h po as 

o _ fviVjdvi ~ 9viVj fvi fviQvk ~ 9vifvk 



5,= 



(vt — vi)'-- (vi - v n )g Vi ' (vi - vi) ■ ■ ■ (vi - v n )g Vi {vi - v k ) ' 

gfvi - fgv, 



9vt 

Conjecture 1 The corresponding (n+A) -dimensional system 

p n+1 n+1 

(Vi)t - ^(Vi) x + ]T 7rf(«0»™ +*i(«i)y» +1 + E Ai«r(«i)»™ = 0, (82) 

^ m=0, m=0 : 

i = 1, . . . , ii, is integrable for any n. 



35 



We note that for n = 3, we can easily reduce system (|52"j) to the form 
analogous to (|8Tt 



"xi - w X2 + uw X3 - wu X3 = 0, v yi - w V2 + vw y3 - wv V3 = 0, 
w z = v Xl - u yi + vw X3 - wv X3 - uw V3 + wu y3 . 

We plan to reduce system (|82l) to a simpler form for arbitrary n, to prove 
that it is integrable, and to study the corresponding GT system in a separate 
paper. Here, we only note that the definition of the GT system as a compatible 
system of form (|28|) can be easily generalized to the case d > 3. Moreover, if 
GT systems in the three-dimensional case are related to fiber bundles with a 
one-dimensional fiber (see SeclHl), then a fiber in the d-dimensional case has the 
dimension d — 2. Some GT systems for d > 3 were found in [15]. We present a 
few new examples of one-field GT systems in the case d = 4. 

Let P(x) = Z2X 2 +Z1X+Z0 and Q(x) be two arbitrary quadratic polynomials, 
J{x) be an arbitrary linear polynomial, and S(x, y) — 2z2xy + z\{x + y) + 2zq. 

Example 8 The formulas 

P(Pi)i.Qi - Qj) P{Pi) 

d ^ u - P( pl )p( pj )( ql - qj ) d * ud > u 

give a one-field GT system with a two-dimensional fiber. The coordinates in the 
fiber are pi and qi, and u is the coordinate on the base. Using a fractional-linear 
transformation ofp, we can reduce the polynomial P to one of the two canonical 
forms P(x) = x or P(x) = 1. The case P(x) — 1 in other coordinates arose 
in when investigating the equation 

Zt X + Z X y + Z XX Zyy ~ Z^y = 0. 

The general position case P{x) = x is possibly related to the most nondegenerate 
integrable equation of this sort. 

Example 9 The formulas 

Q V _ (Pi-Pi) J (Pi)J(Pi) d g J_Mt Q 



dAu ~ pipmPi)^-^) 2 1 3 



give a one-field GT system. 
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Example 10 The formulas 

p%Pj(2-pj) Pi 

didnU = - — — — -^ diudjii 
(ft"*) 2 

give a one-field GT system. 

We do not know whether the GT systems in Examples |H] and [TU] are related 
to some multidimensional integrable dispersionless systems. 

D Classifying integrable chains of hydrodynamic 
type 

We consider integrable quasilinear infinite chains of the form 

U a ,t — 4>a,lUl iX H h (t)a,a+lU a + l. x , a = 1,2,..., 4> a ,a + l ^ 0, (83) 

where <j) a _j — cf> a j(ux, . . . , u a+ i). Two chains are equivalent if they are related 
by a transformation of the form 

tt a -»* a (ui,...,ii a ), - — ^0, a = 1,2,.... (84) 

We consider a chain integrable if it admits hydrodynamic reductions (see |13|- 
EE E3]-[25]). 

Example 11 TTie Benney chain l45^ - \4^i 

Ul tt = U2,x, U2,t = ttl«l,a:-Ht3,X) «at = (a-lK_iUi iS; +U a+ i. r , 

(85) 

is £/ie best-known example of an integrable chain (|83l) . Hydrodynamic reductions 
of the Benney chain were studied in JZ1J/. 



Definition 9 ^4 hydrodynamic (1+1) -dimensional N -component reduction of 
chain (|83l) is semi-Hamiltonian system (|43[) m which the functions Uj (r 1 , . . . , r ), 
j = 1, 2, ... , satisfy system (|83|) /or every solution (|43p. 



Using transformations (|29p . we can set 

F(p,wi, • • . ,«„) =p (86) 

without restricting the generality. 

A calculation analogous to the one at the end of Sec. [5] results in an infinite 
triangular GT system corresponding to the given integrable chain. 
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Definition 10 We call a compatible system of the form 

diPj = f(Pi,Pj>Ui, . . .,U n )diUi, 

diU k = g k {pi,ui,...,u k )diUi, fc = l,2, (87) 

didjUi = h(pi,Pj,ux, . . . ,u n )diUidjUi, 

where i,j = 1,...,N, i + j, a triangular GT system. Here, pi,...,pn and 
Ui, U2, ■ ■ ■ are functions of r 1 , . . . , r N and d{ = d/dr l . 

Substituting m — u^r 1 , . . . , r ), i = 1, 2, . . . , in the chain, calculating the 
derivatives with respect to t and x by virtue of (|43]) , and equating the coefficients 
of r s x to zero, we obtain 

d s u a p s = <j) a ^d s ui H h (/) ai a + id s u a+1 , a = 1 , 2, . . . . 

Using (|87l) and replacing p s with p, we find 

P = <£l,l+01,252i = 02,1 + 02,252+02,353, P93 = 03,1 + 03,252+03,353 + 03,454, 

Solving this system for 52,53, • ■ • , we obtain 

ffi(p) = V»i,o + ipi,iP H 1- V'M-iP 1-1 , 

where are some functions of iti, . . . , ttj. In particular, 

52 = — • (88) 

01,2 01,2 

Example 11a (continuation of Example [TTj) . The triangular GT system 
corresponding to the Benney chain is 



diPj = , didjUi = +y, (89) 

Pi -Pi [Pi-Pj) 

d t u m = (-(m - 2)u m - 2 2u 2P ™- 2 - u lP ™- 3 + p™- 1 ^^ (90) 



Equations (|59"f were first obtained in [T2"] . 

The compatibility conditions didjU a = djdiU a , a = 2, 3, 4, result in a system 
of linear equations for diPj, djPi, and di djU\. Solving it, we obtain 

d iPj = P(Pi,Pi) d iUl , (91) 
Pi - P 3 



Q{pi,p 3 ) 
(pi -Pif 



didjUi = — — % '* 2 diUidjUi, i + j, (92) 



where P and Q are quadratic polynomials in each of the variables pi and pj. 
Their coefficients may in principle depend on u\, U2, U3, and U4. But it is easy 
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to deduce from the compatibility conditions didjPk = djdiPk that P and Q 
depend only on u\ and u%. 
We write (|9"Tj) in the form 

diPj = ( R 3 + {Zip* + z 5 p 3 + z 6 )p t + ZiPj + z^p* + z 7 pj + z s j dim, (93) 

\Pi Pj / 

where R(x) — z 4 x A + z^x 3 + z 2 x 2 + z\x + zq. From the GT-system compatibility 
conditions, we easily obtain 

/ 2z 4 pf p] + z 3 p iPj { Pi +Pj) + z 2 {pj +p]) + Zl ( Pl + Pj ) + 2z \ 

OidjUi = hZg WiUlOjUl. 

V (Pi Pj) J 

(94) 

Using change of variables (|84|) . we can set the coefficient zg equal to zq — zj. The 
coefficients z%{x, y), i = 1, . . . , 8, satisfy a pair of compatible dynamical systems 
in y and x. The first of these systems is 



z 0,y 


— 2zqZ$ 


- z l z 6, 


z i,y 


= 4z z 4 + z ± z 5 - 


- 22226, 




= 3z\Z4. 


— 3Z3Z6; 


Z 3,V 


= 2z 2 z 4 — z 3 z 5 - 


- 4z 4 z 6 , 


z i,y 


= Z3Z4 - 


- 2z 4 z 5 , 


z b,y 


= Z4Z7 — Z 4 Zq — 


z l 


z 6,y 


— Z 4 Z 8 - 


' Z 5 Z 6> 


z T,y 


= 2z\z 4 — 2z^zq 


— Z 5 Zq + Z 4 Zg, 


z s,y 


= 2zqZ 4 


- z\ - z 6 z 7 + z 5 z s . 









The second of these systems looks more complicated. To linearize these systems, 
we reduce the polynomial R to a canonical form, sacrificing normalization (|86l) 
for this. 

It turns out that if the transformation coefficients 

m+b i = 1 N (95) 

Pi-y 

satisfy the conditions 

a« 2 = z i( b + a V0i K 2 = z iH + z 5 b - z 6 a, ipu 2 = z ifp 2 + z 5^ + z 6, 

then the transformation preserves the form of Eqs. (p?3"|) and (|94|) . The polyno- 
mial R changes simply under transformations (|95p: 



R{pi) -»• {Pi - ^) A R 



api + b 
Pi-tp 



We first assume that all the roots of R are distinct. We can then verify that using 
admissible transformations (|95l) . we can set three of the four roots to be 0, 1, 
and 00. The GT-system compatibility conditions then state that the fourth root 
X(m,u 2 ) is independent of u 2 . Using transformations of the form u\ —> q(ui), 



39 



we obtain cither A = u\ or A = const. The GT-system compatibility conditions 
then imply that in the first case, Eqs. (l93|) and (|94|) coincide with the formulas 
in Example^ where P(x) = x{x — 1), while the second case cannot be realized. 

In the case of multiple roots, we can reduce the polynomial R(x) to one of 
the canonical forms R = 0, R = 1, R = x, R = x 2 , or R = x(x — 1). In all 
these cases, Eqs. and (IM1) coincide with the corresponding equations in 
Example [1] 

Below, we consider the general position case R(x) — x[x — l)(x — u\) and 
the most degenerate case R(x) = 0. 

Because we use transformations (|95[) when reducing R to a canonical form, 
the functions gk in system (|87[) change their form from polynomial to rational 
functions with the denominator (p—ip) k ^ 1 . To construct them, we must describe 
all possible fractional-rational functions g 2 . Using a transformation of the form 
U2 = <j{u\, U2) in the general position case, we can bring every such function g 2 
to one of the functions 

u 2 (u 2 - l)(p-ui) . . . 

1- ff2(p) = — 7 77 (the regular extension), 

ui(ui -l)(p- u 2 ) 



2- 92(p) 



1 



p-ui 



3 . B ,0,) = ^>- - A = 1,0, 

P- A 

, . Ml - U 2 M 2 - 1 

4 - 92b) = —7 -tP+ r- 

Mi(Mi — lj Mi — 1 

The regular extension, case 1, is invariant under the discrete automorphisms 
described in Example [2] The other three cases are equivalent, and we can 
consider case 4 as an example. 

We now consider the general position case 1. The next step in the classifica- 
tion is finding fractional-rational GT families with the coefficients depending on 
Mi and u%. Such a general form family is determined by formulas (|50l) and (|5T|) . 
where n = 2 and k = 1. But we have the additional constraint that a zero of 
the denominator must coincide with a zero of the numerator of g 2 , i.e., it must 
be equal m 2 . It is easy to verify that this condition is equivalent to s 2 — and 
h>i,u2 = ^2,u 2 = in system (|T7|) . (14*51) . The last condition implies that hi{u\) 
and /i 2 (wi) are linearly independent solutions of the standard hypergeomctric 
equation 

m(m-1)/i(m)"+[si + S3-(s 3 + S4 + 2si)m]/i(m)'-(-si(si + s 3 + S4 + 1)/i(m) = 0. (96) 
Without restricting the generality, we can choose 

Jo 
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As a result, we obtain 



where 



jp, \ fi(ui,u 2 )p- f%{ui,u 2 ) , , 

F(p,u 1 ,u 2 ) = , (97) 

P- u 2 



u 2 (u 2 - l)hih 3jU2 + ui(«i - l)(ft-i/i 3 .«i - h 3 h[) 
h — 



h = 



u\(u\ — \){hih' 2 — h 2 h'i) 

mu 2 {u 2 - l)hih 3tU2 + u 2 u 1 (ui - l)(/iife3,m ~ fe 3^-i) 
ui(ui — l){hih' 2 — h 2 h'i) 



We note that hxh' 2 - h 2 h[ = const • (w a - l) Sl+S4 ttf +S3 . 

For special values of the parameters, we can solve the hypergeometric equa- 
tion in elementary functions and find F explicitly. For example, we have 

_ (u 2 - u 1 )^+ 1 u s 2 3+1 {u 2 - 
P - u 2 

for s 4 = —2 — Si — s 3 and 

(p ~ 1)(«2 " Wl) S1 + 1 ^ 3 + 1 (ui ~ I)- 1 " 81 



p-u 2 



for S4 = 0. 

Further, we must find the functions g 3 , 34, ... in system (|87j). These func- 
tions are defined up to an arbitrary transformation (1841) . where a = 3,4, ... . 
To fix them more or less rigidly, it is convenient to assume that the functions 
3 3 , 34, . . . are linear in itj, i > 2 (cf. (|90[1 ). In particular, we can choose 

, , {u\- u 2 ){u 2 - l)p 

93{P) = j T77 yj, 

m(ui - i)(p- u 2 y 

^ _ (j - 3)(m - u 2 ){u 2 ~ l)puj _ (ut - u 2 y~ 3 (u 2 - l) 2 p{p - ui)(p - 1)'~ 4 
Ux(u\ - l)(p - u 2 ) 2 ui(ui - iy- 2 (p - u-z) 1 - 1 

^ (i-s- 2)(u 1 - u 2 ) s (u 2 - l) 2 p(p - MiXg - l) s ~ 1 u i _ s 
^ - l) s+1 (p-u 2 ) s + 2 ' 

s— 1 v 7 v ' 

The coefficients <^ j of the corresponding chain (|83[) can be determined from 
the relations 

F (P) = 01,1 + 01,252, ^(P)52 = 02,1 + 02,252 + 02,333, 

(98) 

F{P)93 = 03,1 + 03,232 + 03,353 + 03,4#4, . . . , 
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where F is given by formula (j9"T)l . These relations are equivalent to an infinite 
triangular system of linear algebraic equations. Solving it, we obtain 

u 2 {u 2 - l)(ui - U 2 ) 

f 2 ui - hu\ (99) 
u 2 (ui - u 2 ) ' 

02,3 = f\U 2 - f 2 , 

All the functions <f>i,j, i > 2, are linear combinations of f\ and f 2 with ratio- 
nal coefficients with the same denominator (u 2 — l)(wi — u 2 ) and a numerator 
depending on u\, . . . , iti+i- 

We next consider the most degenerate case R(x) — 0. In this case, the 
GT system is holomorphic on the diagonal. According to conjecture [T] (sec 
Appendix |C|) . precisely the models related to such GT systems admit multidi- 
mensional generalizations. 

We can easily verify that in the case R(x) — 0, the triangular system is 
equivalent to the system 

diPj = 0, d l d ui = 0, diU k = P^ui, k = 2, 3, 

Automorphisms of this system are generated by the transformations 

Pj -+Pj, Ui -> vu { +7i, 

(100) 

Pj^-apj+b, Ui -> a 1 1 u t + (i - l)a l 2 bui- 2 + ■ ■ • + h % X U\, 

where j — 1, . . . , N, i = 1,2,.... The corresponding GT family is F(p) = 
A(u\, u 2 )p + B(u\, u 2 ). The coefficients A(x, y) and B(x, y) are found from the 
scmi-Hamiltonian condition, which is equivalent to the system of differential 
equations 

AByy = Ay By, AB X y = Ay B x , AB ' xx = A X B X , 

(101) 

AAyy = Ay, AA X y — A x Ay , AA XX = A x + A x By — AyB X . 

We can easily solve this system in elementary functions. For each solution of it, 
formulas (j98)) determine the corresponding integrable chain (j83|). 

Formulas (|101|) imply that F can depend on u 2 in two different ways: 

1. F{p,u\,u 2 ) = e Xu2 {a{u\)p+ b{u\)) in the general position case, or 

2. F(p, ui,u 2 ) = a(u{)p + Xu 2 + b(ui). 

In the first case, we in turn have two possibilities: b' ^ and b' = 0. We have 
a = a' , b = k\a, a(x) — c\e^ lX + C2e At22: , where c\c 2 {\ki — [iifi 2 ) = 0, 



flU! - f 2 
u x —u 2 

(U2 - l)(flU 2 - f 2 ) 

(ui - l)(ui - U 2 ) 
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and 

b = Ci, a(x) — C2e^ x + C3, where C2(ciA — C3yu) = 0, 
for 6' = 0. The same variants in case 2 correspondingly give 

a = a 1 , b = ki<j, er(a;) = c\ + C2X + C3e Ma: , where C3(A — c 2 ^) = 0, 

and 

b = ci, a(x) = C2e MK + C3, where C2(A — 03^) = 0. 

It is easy to see that in the general position case, transformation (llOOp reduces 
the function F to the form 

F(p) = e U2+Ul (p - 1) + e" 2 " 111 (p + 1). 

The corresponding integrable chain is 

u*,t = (e" 2+Ul + e U2 - Ul )u fe+1>:c + (e" 2 "" 1 - e" 2+tll H, x , fc = 1,2, . . . . (102) 

This chain has an infinite hierarchy of commuting flows as usual. For instance, 
the next flow is given by 

Uk.r = (e U2+Ul + e U2 - Ul )u k+2 , x + (u 3 - Ul )(e U2+Ul + e u ^)u k+1 , x + 

+ (e" 2+111 - u 3 - 1 ) + e" 2 -" 1 (u 3 - m - l))u fc , s , k = 1, 2, ... . 

In case 2 with C3 = A = and k± = 1, we obtain the chain 

ujt.t = u k +i, x + uiu k , x , k = 1,2,..., (103) 

which is equivalent to universal hierarchy chain |48) . Chain (|103j) is a degener- 
ation of the chain 

u k ,t = Uk+i, x + U2U k}X , k=l,2, (104) 

Following [TS], [JT], we can easily construct (2+l)-dimensional generaliza- 
tions of the (l+l)-dimensional chains described above. Namely, some families 
of the obtained functions F depend linearly on two parameters. We let 71 and 
72 denote these parameters. The corresponding integrable chain 

Uk.t = 7i(^M u i,x + h 4>k,k+iUk+i,x) + 72(^,1"!,* H 1- ^k,k+iUk+i.x) 

is also linear in 71 and 72. We claim that the (2+l)-dimensional chain 

Uk,t = H h 0fe,fe+lUfc+l,a:) + {tpk.lUl^y H h V , k,/fe+l u fc+l, I /) ( l05 ) 

is integrable from the standpoint of the hydrodynamic reduction method. We 
can easily describe these reductions in each actual case. For example, we have 

F(p) = 7ie" 2+M1 (p - 1) + "f 2 e U2 - ul (p + 1) 
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in the general position case. The corresponding (2+l)-dimensional chain (|105[) 
is 

u k ,t = e U2+Ul (u k+hx - u k , x ) + e U2 - Ul (u k +i,y + u k , y ), k = 1, 2, . . . . 
Its hydrodynamic reductions are given by 

r\ = (e" 2+lll ( K - 1) + e u *-^q t ( Pl + l))r< , rj = q^, 

where 

diUk=Pi~ x ui, k=l,2, d i d j u 1 = 0, = 0, 

d _ (gi - oMji + l)(Pj - i)e M1 + frfo ~ l)(Pj + IK" 1 ] (106) 
3 (Pi-ftO(e Ul +fte-»i) 

Gibbons-Tsarev system (| 106[) with two-dimensional fibers (see Appendix [Cj 
deserves a separate study. 

We now consider one of the degenerate cases F — "f\e Ul p + ^{p + u%). The 
corresponding (2+l)-dimensional integrable generalization of chain (| 104|) is 

Uk,t = e Ul Uk+i, x + Uk+i, y + u 2 u k , y , k = 1,2, 

The hydrodynamic reductions are described by 

r\ = [e Ul Pi + qfyi + u 2 )]r*, r\ = t&r*, 

where 

9iWfe = A; = 1,2, . . . , didjUi = 0, diPj = 0, 

We note that the GT system for this (2+l)-dimensional chain is equivalent to 
GT system (TPS]) . 

Conjecture 2 Every integrable chain (|105l) can fee generated by a two-dimensional 
linear family of solutions of system (jl01[) using the construction described above. 

E List of unsolved problems 

Here, we formulate several unsolved problems that from our standpoint seem 
important for better understanding the nature of the GT system. 

Problem 1 To find all solutions of system (|32p of functional equations up to 
an equivalence. 

Problem 2 To describe all two-field extensions of system (|33|) . 

Problem 3 To describe all solutions of Eq. ((3SJ) for GT system (g5J), (gSJ) (of 

Zeas£ /or small n). 

Problem 4 To construct solutions of Eq. (|25[) and i/ie pseudopotential in the 
case of the regular GT-system extension in Example^ 
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